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Abstract 

We point out that the universality of free-fall can be tested by observing surface-gravity changes 
of the Earth. The Earth's inner core is weakly coupled to the rest part of the Earth by mainly 
gravitational forces. If there were a violation of the universality of free-fall, because of their different 
chemical compositions and/or of different mass fractions of binding energies, the inner core and the 
rest part of the Earth would fall at different rates towards the Sun and other sources of gravitational 
fields. The differential acceleration could be observed as surface-gravity effects. By assuming a 
simple Earth model, we discuss the expected surface-gravity effects of violations of the universality 
and experiments to search for such effects by using superconducting gravimeters. It is shown 
that the universality can be tested to a level of 10 -9 using currently operating superconducting 
gravimeters. Some improvements can be expected from combinations of global measurements and 
applications of advanced data analyses. 

PACS numbers: 04.80. Cc, 91.10.Pp 
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I. INTRODUCTION 



The universality of free-fall, stating that every material (point mass) in a gravitational 
field falls at the same rate, is accepted as one of the most fundamental principles in mod- 
ern physics; the Newtonian mechanics and Einstein's general relativity are based on this 
principle. All experiments performed so far support the universality However, theories 
towards the unification of the four forces (strong, weak, electromagnetic and gravitational 
interactions) typically introduce new interactions and predict violations of the principle 
These theories predict new Yukawa-potential type interactions between 
putative charges, such as baryon number, isospin and electrostatic energy density, which 
are functions of proton and neutron numbers in elements. This prediction implies chemical 
composition-dependence of free-fall. Also, many alternative metric theories of gravitation 
predict violations of the universality of free-fall of massive bodies [7]. Therefore, testing 
the universality at a high sensitivity is expected to make a breakthrough in the current 
understanding of physics. Because there are many unknown factors in those theories, no 
precise prediction has been made for the magnitudes of the putative violations; the univer- 
sality should be tested as precisely as possible. Also, because the characters of the putative 
interactions are uncertain, the universality has to be tested for various putative charges, 
using different kinds of test bodies, at different ranges. In addition, to confirm experimental 
results, it should be tested by at least two different experimental methods. Considering 
the importance of testing the universality and the necessity of variety in experimental ap- 
proaches, we present a new method of testing the universality. So far, it has been tested by 



various experiments (e.g. 



and it is verified to a level of 10 13 by laboratory Eotvos-type 



experiments [9] and by Lunar Laser Ranging (LLR) 10]. 

The merits of the proposed method may be as follows: (1) this is the first application 
of the superconducting gravimeters, whose performance has been verified to be sensitive 
and stable in geophysical studies, to tests of the universality, (2) this is the first attempt to 
use Earth's interior as test bodies, and (3) unlike most of the experiments for testing the 
universality, this method has a potential to test the universality of free-fall of gravitational 
self-energy. 

We present the concept of this method in section [Til the theory in section II II I and the 
expected sensitivity in section | 
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II. THE CONCEPT 



In this method, the test bodies are the solid inner core and the rest part of the Earth (the 
liquid outer core, the mantle and the crust). They are in free fall in the gravitational field 
mainly due to the Sun. If there were a violation of the universality of free-fall, there would 
be differential acceleration between the inner core and the rest part of the Earth towards the 
Sun (see section [TTTI for detail). The differential acceleration can be searched for by measuring 
surface-gravity changes. Relative gravity changes can be measured with high resolution by 
using superconducting gravimeters. They have been used to search for translational motions 
of Earth's inner core (the Slichter triplet [111]) (e.g. It is shown that 



they are capable of measuring long period surface-gravity effects. We attempt to apply the 
geophysical tool to test the theory of gravitation. 

In the test bodies, there are mainly two kinds of differences that could result in the 
differential acceleration due to violations of the universality. One is the difference in their 
chemical compositions (the chemical composition-dependent effect), which implies the dif- 
ference in proton and neutron numbers in the elements that compose the test bodies. The 
solid inner core mainly consists of iron and nickel, and its density is approximately pi C ~ 



13000 kg m 3 17|, while the rest part of the Earth is mainly made of lighter elements such 



as silicon oxides, except the liquid outer core (which mainly consists of heavy elements and 
whose density is slightly less than the inner core (p oc ~ 11000 kg m -3 jl3])). The average 
density of the rest part, including the outer core, is approximately 5400 kg m~ 3 . 

The other is the difference in the fractions of gravitational self-energy in the test bodies 
(about —3.7 x 10~ n and — 4.2 x 10~ 10 for the inner core and the rest part, respectively). 
Testing the universality of free-fall of gravitational self-energy can be viewed as a test of the 
strong equivalence principle. 

The chemical composition-dependent effect can be tested using laboratory-size test 
masses. The effect of gravitational self-energy appears significant only in massive test bodies. 

III. EQUATION OF MOTION OF THE INNER CORE 

We assume a simple configuration as shown in figure [U The Earth is revolving around 
the Sun (with angular frequency Ur) in a circular orbit with radius r. The Earth's rotation 
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axis is perpendicular to the ecliptic plane, namely the inclination of the rotation axis by 
about 23.3° is not considered. The x-y plane is set on the Earth's equatorial plane with its 
origin at the center of figure of the rest part of the Earth. 

The Earth's interior can be classified into four parts: the solid inner core, the liquid outer 
core, the mantle and the crust. We assume that the solid inner core is a homogeneous sphere 



(with density p« c and radius rj C ~ 1.22 x 10 6 m 17|) and it is enclosed in the spherical liquid 
outer core (with the average density p oc and outer radius r oc ). We do not consider any 
deformations (such as tidal or rotational deformations). We assume that the mantle and the 
crust are spherical shells with uniform densities, and their centers of figures are coincident; 
their gravitational influence on the inner core is negligible due to Newton's shell theorem. 
The Coriolis acceleration splits oscillations of the inner core into a triplet of periods (the 
Slichter triplet However, for simplicity, we assume that the inner core oscillates only 

along the x-axis (with angular frequency u ) and we ignore the Coriolis acceleration. We do 
not consider any electromagnetic effects. 

The effective viscosity of the outer core r\ is not well determined and estimates from 
various methods vary from ~ 10~ 3 Pa s to ~ 10 12 Pa s 18]. The Reynolds number (= 



p oc vr ic /ri, where v is the velocity of the inner core) is less than unity when r\ is larger than 
~ 10 7 Pa s and the amplitude of oscillations X (= v/ojq) is nominally 1 m. We assume that 
the friction between the inner core and the outer core is proportional to the velocity of the 
inner core. 

Under these assumptions, the x-component of the equation of motion of the inner core 
(mass rrii C = 4 irrf c pi C /3 ~ 1.0 x 10 23 kg) can be written approximately in a form of so-called 
damped forced oscillation: 



x ~ — uJqX — 2kx + AacosujRt (1) 



where 



w 2 ~ *G Pie Poc p oc « 4.7 x 10-V 2 « {2tt(2.5 h)" 1 } 2 (2) 
3 Pic 

2k = ~ 2.3 x 1(T 16 77 s- 1 (3) 

The stiffness Uq is mainly due to the gravitational pull by the outer core and, as seen in 
equation ([2]), determined by the density difference in the core with the gravitational constant, 
G = 6.67 x 10~ n N m 2 kg -2 . The stiffness due to Sun's tidal force that acts to enlarge any 




FIG. 1: A schematic cross section of the assumed configuration (not drawn to scale). The Earth is 
fixed and the Sun goes around the Earth in the circular orbit (radius r) at the angular frequency 
ujr ~ 27r(24 h) . The inclination of Earth's rotation axis to ecliptic is assumed to be 0° for 
simplicity. The x-y plane is set on the Earth's equatorial plane. The Earth's inner core oscillates 
along the x-axis. The gravitational influence from Earth's mantle and crust (not drawn) is negligible 
because of Newton's shell theorem (see text). 

displacements of the inner core from the center of the Earth twice per its revolution is about 
8 orders of magnitude smaller than the gravitational stiffness and ignored. The stiffness due 
to the Moon's tidal force is about twice larger than the one due to the Sun's tidal force and 
ignored. The stiffness due to the centrifugal force caused by Earth's revolution (« cj^) is 
about two orders of magnitude smaller than the gravitational stiffness and ignored. 

Aa in equation flTJ is the magnitude of the putative differential acceleration due to a 
violation of the universality of free-fall in the Sun's gravitational field. The signals of the 
violation are expected to have the frequency of once per revolution of the Earth: ujr ~ 2-7r(24 
h) _1 m 7.3 x 10 -5 s _1 . At this frequency, the major obstacle would be the tides (the body 
tides and the ocean tides) which have the same frequency as the violation signals. The tidal 
effects and other known effects can be largely removed from the gravity data by applying 
appropriate models and data analysis methods (e.g. [3]). 
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IV. EXPECTED SENSITIVITY 



Assuming that k < ojq (or 77 < 5.9 x 10 12 Pa s), the solution of equation p]) can be 
written by a sum of two terms that express damped oscillations and forced oscillations. The 
damped oscillations decay exponentially with time t (e~ kt ; k~ x is about 24 minutes for rj = 
5.9 x 10 12 Pa s or when k = ujq), while the forced oscillations remain: 

x = Asin(u B t - 5) (4) 

where 

A = Aa (5) 

+ x 2kujR (a\ 

tan6 = —2 —2 ( 6 ) 
uj - u R 

4 

When the damping coefficient is sufficiently small (i.e. 4k 2 C ^- ~ 4 x 10~ 5 s~ 2 or rj <C 

LO R 

2.8 x 10 13 Pa s), we obtain 

and tan 5 ~ 0. When the inner core is shifted from the center of the Earth by A max , along 
the x-axis, this displacement will change the gravitational field at a surface point on the 
x-axis by 

. 2Gm ic Aa Aa 
r\ u 2 10 

where r e = 6.4 x 10 6 m is the average radius of the Earth. 



The noise leve 



of superconducting gravimeters at quiet sites is a few 10 m s at the 



signal frequency Assuming a somewhat better sensitivity of 10 12 m s 2 at the signal 
frequency, we could estimate the universality of free-fall to 

Aa 



-3 ™ c-2 i c 



1.6 x 10~ 9 (9) 



where a sun ~ 5.9 x 10 m s is the Sun's gravitational acceleration that acts on the Earth. 
Because of the inclination of the Earth's rotation axis, the maximum violation signals can 
be expected at observatories located on the equator in Spring and Autumnal equinox points, 
and on Tropic of Cancer or Capricorn in Summer and Winter solstices. 

Equations ([7]) and ([8]) indicate that the relative displacement of the inner core can be 
measured to ~ 20 /mi. For this amplitude, the Reynolds number is unity when r] ~ 200 Pa 



s. As a result, the range of the effective viscosity, discussed above, is 200 Pa s < rj < 6 x 
fO 12 Pa s. This range is consistent with the estimates summarized in 

V. DISCUSSION 

This work is based on the simple Earth model and configuration. Elaborate modellings 
have to be applied for more accurate analyses of expected violation signals and of the ex- 
pected sensitivity Especially, the value of which affects the estimates of the expected 
signals and sensitivity, is highly model dependent. The latest theoretical studies with elab- 
orate Earth models predict smaller values of ivq than the value we have used (equation (T5])), 



between ~ 2n(6 h) 1 and ~ 27r(4 h) 1 [20l . |21| . Because the expected sensitivity is pro- 

X 2 , 



portional to ~ cu 2 , the expected sensitivity could be better by ~ 3 to 6 times by applying 



elaborate Earth models. However, we have assumed the rather optimistic sensitivity of 10~ 12 
m s -2 . Therefore, the expected sensitivity would remain to be on the order of 10 -9 . 

The current limits on violations of the universality of free-fall are on the order of 10~ 13 . 
In order to get comparable results, a sensitivity better than a few 10~ 15 m s~ 2 (0.1 picogal) 
at the signal frequency is required. Presently, superconducting gravimeters are the most 
sensitive instruments for measurements in the frequency range. Though it seems difficult to 
achieve this sensitivity, there are several possibilities to improve the sensitivity. One way 
may be applications of more sophisticated data analyses than the usual Fourier analyses, 
as discussed in [22] to search for the Slichter triplet. Another way may be to carry out 
coincidence measurements with two superconducting gravimeters located ideally opposite 
sides of the Earth near the equator. If there were a violation towards the Sun, the expected 
magnitude of the violation signal at the two superconducting gravimeters is the same but 
the sign should be opposite. By combining such coincidence signals, we could double the 
magnitude of the expected signals and the sensitivity would be improved by a factor of two. 
Candidate observatories for such coincidence measurements may be the one at Hsinchu in 
Taiwan (25°N 121°E), where we are currently setting up two new superconducting gravime- 
ters, and the one at Concepcion in Chile (37°S 73° W). If the noise level of data from the 
observatories were high, it might be better to use data from low noise sites considering the 
degrees of signal compensation depending on the latitude and longitude of the sites. Such 
global observations would be possible through the Global Geodynamics Project network 
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(GGP |23j). Further studies are necessary to figure out the optimal schemes for global 
observations and noise reduction. Currently, we are checking the performance of the new 
superconducting gravimeters, installed in Hsinchu Taiwan in March 2006. We plan to carry 
out a preliminary test of the universality in the near future. 

As described in section [Til the difference in the mass fraction due to gravitational self- 
energy is a grav ~ 4 x 10~ 10 . Therefore, from equation ([9]) and the above discussion, the 
geophysical test with current superconducting gravimeters would be only sensitive to the 
chemical composition-dependent effect, but not sensitive enough to test the universality of 
gravitational self-energy. If the sensitivity were improved to ~ 0.1 picogal, we could test the 
univers ality of free-fall of gravitational self-energy to the same level as the LLR experiment 



id, 



24 



25]. 



VI. CONCLUSIONS 



We have considered a new method of testing the universality of free-fall. In this method, 
differential acceleration between the Earth's inner core and the rest part of the Earth is 
to be searched for by measuring surface-gravity effects. Based on a simple model, we have 
shown that the universality would be tested to a level of 10~ 9 with current sensitive super- 
conducting gravimeters. Some improvements can be expected from combinations of global 
measurements and developments of methods of data analysis. We plan to carry out a pre- 
liminary test of the universality using superconducting gravimeters in Hsinchu Taiwan. To 
get a comparable result with the LLR experiment, the sensitivity has to be improved by 
about four orders of magnitude at the signal frequency (once per day). A breakthrough in 
developments of gravity measurements is necessary to achieve this sensitivity. 
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